Time drift of cosmological redshifts as a test of the Copernican principle 
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We present the time drift of the cosmological redshift in a general spherically symmetric space- 
time. We demonstrate that its observation would allow us to test the Copernican principle and 
so determine if our universe is radially inhomogeneous, an important issue in our understanding of 
dark energy. In particular, when combined with distance data, this extra observable allows one to 
fully reconstruct the geometry of a spacetime describing a spherically symmetric under-dense region 
around us, purely from background observations. 

PACS numbers: 98.80.-k, 98.80.Es, 95.36.+X 



I. INTRODUCTION 

Cosmological data is usually interpreted under the as- 
sumption that the universe is spatially homogeneous and 
isotropic. This is justified by the Copernican principle, 
stating that we are not located at a favoured position in 
space. Combined with the observed isotropy, this leads 
to a Robertson- Walker (RW) geometry [J, at least on 
the scale of the observable universe. 

This implies that the spacetime metric reduces to 
a single function of the cosmic time, the scale factor 
a(t). This function can be Taylor expanded as a(t) = 
a + H (t — t ) — ^q HQ(t — to) 2 + . . . where H is the Hub- 
ble parameter and go the deceleration parameter. Low 
redshift observations Q combined with the assumption 
of almost flatness of the spatial sections, justified mainly 
by the cosmic microwave background data [H, 0] , lead to 
the conclusion that qo < 0: the expansion is accelerating. 
This conclusion involves no hypothesis about the theory 
of gravity or the matter content of the universe ||, as 
long as the Copernican principle holds. This has stim- 
ulated a growing interest in possible explanations [H, Q , 
ranging from new matter fields dominating the dynamics 
at late times to modifications of general relativity. 

While many tests of general relativity on astrophysi- 
cal scales have been designed [7], the verification of the 
Copernican principle has attracted little attention, de- 
spite the fact that relaxing this assumption may be the 
most conservative way, from a theoretical perspective, of 
explaining the recent dynamics of the universe without 
introducing new physical degrees of freedom [8| . 

This possibility that we may be living close to the 
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center (because isotropy around us seems well estab- 
lished observationally) of a large under-dense region 
has attracted considerable interest. In particular, the 
low redshift (background) observations such as the 
magnitude-redshift relation can be matched @ by a non- 
homogeneous spacetime of the Lemaitre-Tolman-Bondi 
(LTB) family (that is, spherically symmetric solution of 
Einstein equations sourced by pressureless matter and no 
cosmological constant). Unfortunately, this simple exten- 
sion of the RW universes depends on two free functions 
(see below for details) so that the reconstruction is under- 
determined and one must fix one function by hand. Thus, 
one needs at least one extra independent observation to 
reconstruct the geometry of an LTB universe. A limita- 
tion to this reconstruction arises because most data lie 
on our past light cone. This takes us back to the obser- 
vational cosmology program [l(| and the question of 
how to extract as much information as possible about our 
spacetime from cosmological data alone. Among many 
results, it was demonstrated [l2j that the two free func- 
tions of a LTB spacetime can be reconstructed from the 
angular distance and number counts, even though evolu- 
tion effects make it impossible to be conclusive |13| . 

Recently, two new ideas were proposed. First, it was 
realised [14[ that the distortion of the Planck spectrum 
of the CMB allows one to test the Copernican principle. 
Second, a consistency relation between distances on the 
null cone and Hubble rate measurements in RW universes 
was derived [15j . based on the fact that the curvature is 
constant; this also serves as an observational test of the 
Copernican principle. 

In this letter, we reconsider the time drift of cosmolog- 
ical redshift in spacetimes with less symmetries than the 
RW universe and we demonstrate how, when combined 
with distance data, it can be used to test the Copernican 
principle, mainly because observing the thickening of our 
past light cone brings new information. As pointed out 
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by Sandage and then McVittie [16[, one should expect 
to observe such a time drift in any expanding spacetime. 
This may lead to a better understanding of the physical 
origin of the recent acceleration [l?], fl8l. or to tests of the 
variation of fundamental constants [191 ]. This measure- 
ment, while challenging, may be achieved with Extremely 
Large Telescopes (ELT) and in particular, it is one of the 
main science drivers in design of the COsmic Dynam- 
ics Experiment (CODEX) spectrograph [2(J. Our result 
may strengthen the scientific case for this project. 

We start by introducing observational coordinates, 
which allow us to derive the general expression for the 
time drift in a spherically symmetric, but not necessarily 
spatially homogeneous, universe [see Eq. ((6])]. We show 
that observation of both the luminosity distance and the 
redshift drift allows one to probe the Copernican princi- 
ple at low redshifts, when "dark energy" dominates [see 
the consistency relation ©]. We demonstrate that this 
expression may be used with distance data to fully re- 
construct the geometry of a LTB spacetime [see Eq. ©] . 



II. i IN A SPHERICALLY SYMMETRIC 
UNIVERSE 

Observational coordinates. We consider a spher- 
ically symmetric spacetime in observational coordinates 
{w, y, t9, tp}, where w labels the past light-cones of events 
along the wordline C of the observer, assumed to lie at 
the center so that w is constant on each past light cone, 
with u a d a w > 0, u a being the 4- velocity of the cosmic 
fluid, u a u a = — 1. y is a comoving radial distance co- 
ordinate specified down the past light cone of an event 
O on C. Many choices are possible such as the affinc 
parameter down the null geodesies from O, the area dis- 
tance, or the redshift; whatever choice we assume on the 
past light cone of 0, it is specified on other past light 
cones through being comoving with the cosmic fluid, i.e. 
y. a u a = 0. (t9, ip) are angular coordinates based at C 
and propagated parallelly along the past light cone. The 
metric in these observational coordinates is 



ds 2 = -A 2 (w,y)dw 2 + 2A(w,y)B(w,y)dydi 
+C 2 ( Wl y)dn 2 , 



(1) 



which is clearly spherically symmetric around the world- 
line C defined by y — 0. The requirement that the 
2-spheres {w,y} = const, behave regularly around C 
when y — > implies that A(w,y) — > A(w, 0) ^ 0, 
B(w, y) B(w, 0) ^ and C(w, y) = B(w, 0)y + 0{y 2 ). 

There remain two coordinate freedoms in possible 
rescalings of w and y. However, once specified on C, 
w is determined on the other worldlines by the condition 
that {w = const.} are past light-cones of events on C. 
This allows us to arbitrarily choose A(w,0). Also, once 
specified on one past light cone, y is determined on all 
the others because it is a coordinate comoving with the 
fluid. This allows us to choose B(wo, y) for a given value 
of w — Wq. 



On each past light cone, the cross-sectional area of a 
source is related to the solid angle dfl 2 under which it is 
observed by an observer onC a,tw = wq by C 2 (wo, y)dfl 2 . 
This implies that C is the angular distance, Da, i.e. 
DA(y) = C(wQ,y). The distance duality relation [2l[ 
then implies that the luminosity distance is given by 
Dl{v) = (1 + z) 2 Da- The redshift is given by 
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(U a k a ) 



emission 



(u a k a ) 



observer 



A(w ,0) 



(2) 



where the matter velocity and photon wave-vector are 
given by u a = and k a = (AB)~ 1 6y respectively. 

We deduce that the isotropic expansion rate, defined by 
3H = V a u a , is given by 



d w B{w, y) ^ d w C{w,y) 
B(w,y) C(w,y) 



(3) 



For the central observer, who sees the universe isotropic, 
H is simply the Hubble expansion rate. At small red- 
shifts, H(w,y) = B (Zfi)A(w,o) +°(y)' so the Hubble con- 
stant is H = d w B(wo,0)/B(wo,0)A(wo,0). 

In the particular case of a dust dominated universe, the 
acceleration and vorticity vanish and the fluid 4- velocity 
can be expressed as the gradient of the proper time along 
the matter worldlines: u a = —d a t. Since we also have 
u a = —Ad a w + Bd a y we deduce that dt = Adw — Bdy 
so that A = d w t and B = d y t. The surfaces of simul- 
taneity are thus given by Adw = Bdy and we have the 
integrability condition d y A + d w B = 0. 

The covariant derivative of u a is therefore of the form 
V a u b = H(g ab + u a u b ) + a ab , where the shear a ab is 
symmetric, traceless, and satisfies u a a ab = 0. The 
scalar shear a 2 — cr ab a ab /2 is consequently the only non- 
vanishing kinematical variable and is given by 



a(w,y) = 



_2_1_ ( d w B d w C 
V3A C~ 



(4) 



where an arbitrary sign has been chosen. The regularity 
conditions imply cr(w, 0) = 0, which is expected since the 
expansion is observed to be isotropic about the central 
worldline. 

Expression of the redshift drift. From the ex- 
pression ([2]), it is straightforward to deduce that i = 



^(w ,y) is given by 



z(wo,y) = (1 + 2) 



d w A(w o ,0) d w A(w ,y) 



A(w a ,0) A(w ,y) 



.(5) 



Now, we can choose w such that A(wo,0) = 1. Then, 
on our past light cone, we can choose y such that 
d w In B(wq, y) — d w h\A(wo,y). Note however that such 
a choice is not possible for all w. It follows that 



z{wo,y) = (1 + z ) H o - H(w ,y) 



1 



V3 



a(w ,y) . (6) 



This is the general expression for the time drift of the 
redshift as it would be measured by an observer at the 
center of a spherically symmetric universe. 
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Robertson- Walker case. Since a — for a RW 
spacetime, Eq. Q reduces to the standard Sandage- 
McVittie formula. Let us consider the RW metric in 
conformal coordinates 



d s 2 =a 2 (77) [-d77 2 



d x 2 + f K (x)dn 2 ] 



(7) 



where x is the radial comoving coordinate and fx = 
(sin x, X) sinh x) according to the curvature of the spatial 
sections. Now, with w = r/ + x an d y = X> this leads to 
the form (fl]) for the metric (J7J) with A = B = a(w — y) 
and C = a(w — y)fi((y)- H is thus constant on each 
constant time hypersurface and a = everywhere. Since 
d w A\ y=const . = d n a, Eq. © gives z = (1 + z)i?o - #(2) 
where we have shifted to cosmic time (dt — adrj), using 
1 + z — ao/a. 

We understand from this exercise why the observa- 
tional coordinates are adapted to the computation of z in 
an arbitrary spacetime. They are just a generalisation to 
less symmetric spacetimes of the conformal coordinates, 
where the expression is easily obtained in the RW case. 

We also conclude that since the 3 functions (A, B, C) 
are expressed in terms of a unique function, all back- 
ground observations depend on some function of H , z 
being no exception. 

Consistency relation. It follows that, in a RW uni- 
verse, we can determine a consistency relation between 
several observables. From the metric (O and the rela- 
tion for x( z ) that follows, one deduces that H~ 1 (z) = 

D'(z) [1 + n K0 H$D 2 (z)] ' , where a prime stands for 
d z and D(z) = Dl(z)/(1 + z). This relation is the basis 
of the test of the RW structure, as recently proposed in 
Ref. [l5j], arguing there that knowledge of H(z) at dif- 
ferent redshifts, from e.g. baryon acoustic oscillations or 
differential age estimates of passively evolving galaxies, 
could then be used to check this yields the same value of 
VLko- Here, we argue that it can be implemented using 
z(z) as an observational input. Defining 

Cop[D L (z),z(z),z] = l + n K0 H 2 D(z) 2 

- [H (l + z)-z(z)] 2 [D'(z)f(8) 

we must have Cop[-Dl(z), z(z), z] = whatever the mat- 
ter content of the universe and the field equations, since 
it derives from a purely kinematical relation that does 
not rely on the dynamics (i.e. the Friedmann equations). 
This is a consistency relation between independent ob- 
servables that holds in any Robertson- Walker spacetime. 

Spherically symmetric spacetimes. Writing the 
LTB metric in observational coordinates requires the so- 
lution of the null geodesic equation, which is in general 
possible only numerically. Consider an LTB spacetime 
with metric 

ds 2 = -dt 2 + S 2 (r, t)dr 2 + R 2 (r, t)dtf 



where S(r,t) = R'/y/l + 2E{r) and R 2 = 
2M(r)/R(r,t) + 2E(r), using a dot and prime to 
refer to derivatives with respect to t and r. The 



Einstein equations can be solved parametrically as 
{R(r,r)),t(r),r)} = {^S'fa), T (r) + pfjf^ 

where $ is defined by $ (rf) — (sinh 77 — 77, ?y 3 /6, 77 — sin 77) , 
and £(r) — (2E,2,—2E) according to whether E is 
positive, null or negative. 

This solution depends on 3 arbitrary functions of 
r only, E(r), M(r) and To(r). Their choice deter- 
mines the model completely. For instance (E,M,Tq) = 
(— Kqt 2 , Mor 3 , 0) corresponds to a RW universe. One 
can further use the freedom in the radial coordinate to 
fix one of the three functions at will so that one effectively 
has only 2 arbitrary independent functions. Assume we 
fix M(r). We want to determine {E(r),T (r)} to repro- 
duce some observables on our past light cone. This can 
be represented paramctricaly as {r(z), E(z), Tq(z)}. 

Let us sketch the reconstruction and use r as the inte- 
gration coordinate, instead of z. Our past light cone 
is defined as t = t(r) and we set lZ(r) — R[t(r),r]. 
The time derivative of R is given by R[t(r),r] = 1Z\ = 
^2M Q r 3 /TZ(r) + 2E(r). Then we get R'[i(r),r] = 
TZ 2 {r) = -[K(r) - 3(f(r) - T (r))Wi(r)/2]B / /£? - 
TZi (r)Tg(r) + TZ(r)/r. Finally, more algebra leads 
to R'[t(r),r] = TZ 3 (r) = [TZ^r) - 3M r 3 (t(r) - 
T {r))/n 2 (r)]E'{r)j2E{r) + M r 3 T^(r) /TZ 2 + fti(r)/r. 
Thus, R, R' and R! evaluated on the light cone are just 
functions of TZ(r), E(r), Tq(t) and their first derivatives. 
Now, the null geodesic equation gives that 



dt 
dr 



^2 (r) 



y/l + 2E(r)- 



dz 
dr 



l + z 



n 3 (r), 



and 



d^. 
dr 



1 - 



U 2 (r) 



These are 3 first order differential equations relating 5 
functions 7Z(r), t(r), z(r) E(r) and T (r). To reconstruct 
the free functions we thus need 2 observational relations. 
IZ(z) = Da(z) is the obvious choice. Then, from z(z), 
we have the new relation 





(n 3 


^i\ 2 






TZj 



z - (1 + z)H + 



^3 
TZ 2 



£1 
J 1Z 



-1 2 



(9) 



It was shown Q that the observed Dl(z) can be repro- 
duced from the function To(r) assuming E — 0, or from 
the function E(r) assuming To = 0. Here, we have shown 
that {E(r),T (r)} can be completely reconstructed from 
the data without assumptions. 



III. DISCUSSION 

In this letter, we have shown that observation of both 
the luminosity distance and time drift of the redshift as a 
function of z allows one to construct a test of the Coper- 
nican principle. We have derived the general expression 
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for z(z) in a spherically symmetric spacetime [see Eqs. ([5]) 
and |J5))]. This extends the standard computation which 
was restricted to RW spacetimes and was extended to 
almost RW spacetimes only recently [Hf. 

As a byproduct, this also allowed us to derive the con- 
sistency relation between the observed Dl{z) and 
z(z), thereby extending the result of Ref. [l5| to an al- 
ternate observable. That is, while H{z) characterises 
the local isotropic expansion rate, z(z) gives access to 
the expansion between us and a source. In RW models 
these are trivially related but in general the shear enters 
these observables differently thereby presenting a test of 
the Copernican principle using only background observa- 
tions. We have shown that we can extract the shear as 
a function of z and demonstrated that it allows one to 
close the reconstruction problem for a LTB spacetime. 

Dl(z) can be measured from the observation of 
Type la supernovae, particularly with actual projects 
such as JDEM, up to redshifts of a few. z(z) has a 
typical amplitude of order 5z ~ —5 x 10~ 10 on a time 
scale of 5t = 10 yr, for a source at redshift z = 4. This 
measurement is challenging, and impossible with present- 
day facilities. However, it was recently revisited [23| in 
the context of ELT, arguing they could measure velocity 



shifts of order Sv ~ 1 — 10 cm/s over a 10 year period 
from the observation of the Lyman-a forest. It is one 
of the science drivers in design of the CODEX spectro- 
graph (23| for the future European ELT. The study of the 
precision to which we can check the Copernican principle 
with these two data sets is beyond the scope of this let- 
ter. Indeed, many effects, such as proper motion of the 
sources, local gravitational potential, or acceleration of 
the Sun may contribute to the time drift of the redshift. 
It was shown [22jj, however, that these contributions can 
be brought to a 0.1% level so that the cosmological red- 
shift is actually measured. 

Future high precision data may thus allow a test of 
the Copernican principle, even though observations are 
localized on our past light cone. While important in its 
own right for understanding the foundations of our cos- 
mological model, it is also critical for our understanding 
of the acceleration of the universe - it will permit us to 
be confident that any such acceleration is not simply a 
misinterpretation of the data because of incorrectly as- 
suming the geometry of our universe at low redshift. 
Acknowledgements: We thank B. Bassett, F. 
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